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I recently proposed a method of bosonization valid for systems of an even 
number of fermions whose partition function is dominated at low energy by 
bosonic composites. This method respects all symmetries, in particular fermion 
number conservation. I extend it to treat odd systems and excitations involving 
unpaired fermions. 
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. 1. Introduction 

^ ■ There are many ways to introduce composite bosons as elementaty degrees 

of freedom in fermion systems, but most of them violate fermion number 
' conservation, which is an important symmetry in atomic nuclei and other 

' finite systems.^ I recently developed a method to bosonize even systems 

which respects all symmetries, including fermion number conservation.^ 
■ This method was applied to relativistic field theories at zero'^ and nonvan- 

. ishing* fermion number. I present here an application of this technique to 

I nonrelativistic systems. In Section 2 I make an outline of the formalism, 

giving in Section 3 its physical interpretation. In Section 4 I report the 
effective action of composite bosons and finally in Section 5 the effective 
action of composite bosons plus unpaired fermions. 



2. Outline of the approach 

Consider the partition function of a fermion system 



Tt^ i cxp 



-^{Hp - Hf fip) 



(1) 
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where Tr''^ is the; trace over the Fock space, T the temperature, fip the 
chemical potential and hp the fermion number operator. A sector of np 
fermions can be selected by the constraint 

T^\nZ = nF. (2) 

A functional form of Z can be found by performing the trace over coherent 
states 

|c)=exp^^c]„c„j |0), (3) 

where m are the fermion quantum numbers, their canonical creation 
operators and Cm Grassmann variables. Coherent states satisfy the basic or 
defining equations 

Cm\c) = Cm\c) (4) 

where Cm are canonical destruction operators. In terms of these states I can 
write the identity in the fermion Fock space in the form 

2:= / \{dcl,dcm{c\c)-'\c){c\. (5) 

Using the above resolution of the identity the partition function can be 
written 

Z = Tr^{2:exp -T(i?F - Mf n^)] (6) 

where I introduced Euclidean time spacing r and number of temporal sites 
Lq by setting 

T=-L. (7) 
The trace can be evaluated exactly at the relevant order in r with the result 

Z = j [dc*dc] exp [-Sf{c\c)] . (8) 
Sp is the fermion action in Euclidean time 

Sf = tY^ {clVtCt-l + HF{c;,Ct-l) - I^F C*tCt_i} , (9) 

t 

with the discrete time derivative 

^tf=-{ft+i-ft). (10) 

T 
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In a superconducting system one has to perform a Bogoliubov transforma- 
tion which introduces quasiparticles creation-annihilation operators 

a = Ri{c-Bci) 

a'' = {c'' -cB'')R^ , (11) 

where 

R= + . (12) 

They satisfy canonical commutation relations for an arbitrary matrix B. To 
write down the Bogoliubov transformation of coherent states I introduce 
some definitions. First I define the creation operator of a bosonic composite 
with quantum number K 

^ mi,m2 

Bk is the structure function of the composite and O its index of nilpotency 
(assumed to be independent of K), defined as the largest integer such that 

{BkY^^O. (14) 
Next I expand the matrix B on the matrices Bk 

Lastly I define new coherent states 

\a,b) = exp ( -^amttm 1 \b) = exp(-a • )|6) (16) 

m / 



where 

|&) =exp(6 -5^)10) (17) 

is the quasiparticle vacuum. The Bogoliubov transformation changes the 
coherent states |c) according to 

5|c) = (6|&)-5|c,&) (18) 

because^ 

{c,b\c,b) =cxp(c*c){b\b) . (19) 

I can now write a new expression of the identity in the fermion Fock space 
in terms of the transformed coherent states 

I{b*,b)= [[dc*dc]{c,b\c,b)-'^\c,b){c,b\ (20) 
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and use it to get a new expression of the partition function 



Z = Tr' 



'{l{b*,b) 



exp 



}' 



(21) 



Evaluation of the trace yields an action of the system strictly equivalent to 
the original one, but whose individual terms do not respect several symme- 
tries, in particular fermion number conservation. 

There arc many cases in which all symmetries can be enforced term by 
term proceeding in a slightly different way. I introduce at each Euclidean 
time an independent Bogoliubov transformation <St. I do this by letting the 
expansion coefficients bx to depend on time, bx bK,t while keeping the 
basis matrices Bk fixed. Since nothing depends on the B and therefore on 
the expansion coefficients bK,t I can integrate over the latter ones in the 
partition function with an arbitrary measure dn{b* ,b) 



j dii{b*,b) Tr^ I lll{b*, bt) exp [-t{Hf - np) 



(22) 



Requiring the variables b*,b to transform in the proper way all symmetries 
are restored. 

Performing the trace, which can be done exactly^ at the relevant order 
in r, I get a new functional form of the transfer matrix 

Z = J dn{b*,b) H [da*^^,dam,t] exp [-S{b*,b, a*, a, B^B)] . (23) 



3. Physical interpretation 

The last expression of the partition function contains, in addition to the 
Grassmann variables a*, a, the bosonic variables 6*, 6. Is it possible to as- 
sociate these variables to physical bosons? 

To answer this question let us analyze the nature of the states \a,b). 
They are constructed in terms of quasiparticles and composite boson cre- 
ation operators a^,B\^. I call these states coherent because they are ob- 
tained by a unitary transformation from truly coherent states of fermions 
and they share with coherent states of elementary bosons the property of a 
fixed phase relation among components with different number of compos- 
ites. But the basic property of coherent states cannot be fulfilled 



BK\b) ^ bK\b) 



(24) 
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This is a consequence of the composites commutation relations, which are 
not canonical 

1 



Bj,B 



Tr (BjBI 



(25) 



In states with a number of composites n << fi, the above equations can 
approximate the canonical ones by an appropriate normalization of the 
structure functions provided they are sufficiently smooth. Indeed in such 
a case the last term is of order n/Q. But in states with n ^ this is 
impossible even with an absolute freedom about the form of the structure 
functions (which are instead determined by the dynamics). The best we can 
do^ is to satisfy them for states with n + k composites, for fixed n ^ Vl and 
|fc| << 51. The above considerations can be made more precise, also taking 
quasiparticles into account, by showing that under the above conditions 
and for a number of quasiparticles much smaller than O 



n 

K,t 



2ni 



I{b\b)\ 



{aljK..{alJ-{bKj''<^...{bKj'<^) 



..6r 



(26) 



Only for such states can the holomorphic variables b*j^,bK be interpreted as 
the fi,elds of bosons with quantum numbers "K" in the presence of fermions 
with quantum numbers "m" associated to the Grassmann variables am- Re- 
stricting myself to such states I can assume the measure 



dn{b*,b) 



n 

K,t 



db*j^fdbK,t 



2ni 



(27) 



The property 



d|6) = (28) 

ensures that there is no double counting: There are no quasiparticles in the 
bosonic composites. 

What to do with states which do not satisfy the above requirements? 
We should integrate them out (which would require a different choice of 

the measure) leaving only states which have a physical interpretation. But 
in some cases we can reasonably assume that we can ignore them, in the 
spirit of a variational calculation. 



4. The effective action of composite bosons 

In a system of fermions whose low energy excitations are dominated by 
fermion composites I can restrict the trace to these composites neglecting 
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the quasiparticles. The restricted partition function can be written 



Zc = Tr^ <Vcexp 



--{Hp- llF flp) 



(29) 



where Vc is a projection operator in the subspace of the composites for 
which I assume the approximate form 

dh* db 



Vc 



2TTi 



m-'\b){b\. 



(30) 



To proceed further I write the fcrmion Hamiltonian in the form 

K 

The one-body term includes the single-particle energy with matrix e, the 
fermion chemical potential /i p and any single-particle interaction with ex- 
ternal fields included in the matrix Ai 

hQ = e- fj.F+M. (32) 

The matrices Fk are the form factors of the potential, normalized according 
to 

tr(4^FKj = 2fi 

6kiK2- (33) 

Now the trace in the partition function can be evaluated ^ yielding its 
functional form 

'db*db] 



Zc 



2m 



exp(-Sc(6*,6,Bt,B)) 



(34) 



in which 



Sc{b*,b,B\B)=T^\iv[-\n[\ + T R" B^VtB] 
+2ii" B^hB -J29K [{R" - 1) Fj,FK 

K 

+{R" B'^f]^) itr(i?" FkB) - R" B'^f]^R' FkB | . (35) 

In this expression "tr" is the trace over fermion quantum numbers, 

h = ho-J29KFj,FK (36) 

K 

and the variables b*,b must be understood at times t,t — 1 respectively, 
unless otherwise specified. So, for instance 

-1 



i?" ^ (l + BlBt-i 



(37) 
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The latter function should be distinguished from the function R introduced 
in Eq.(12) 

R={l + BB^y' ^ (l + Bt-iBiy' = Rt,t_^. (38) 

The derivation of the above functional form of the partition function is 

based only on the physical assumption of boson dominance and the ap- 
proximation adopted for Vc ■ The action Sc is a functional of the structure 
matrices Bk which can he determined by a variational calculation. 

In many-body physics it is often preferred the Hamiltonian formalism. 
The Hamiltonian of the effective bosons, Hb, cannot be read directly from 
the effective action, because Sc does not have the form of an action of 
elementary bosons. Indeed it contains anomalous time derivative terms, 
anomalous couplings of the chemical potential and nonpolynomial interac- 
tions, which are all features of compositeness. Therefore it has been nec- 
essary to devise an appropriate procedure to derive Hb, which is given^ 
in terms of canonical boson operators b^,b, (not to be confused with the 
composite operators W ,B), so that 



-^{Hb- I^B ub) 



(39) 



Tr is the trace on the boson Fock space, /Ub is the boson chemical potential 
and fiB the boson number operator. 

Hb has a closed form but, for a practical use, it is necessary to perform 
an expansion in inverse powers of the index of nilpotency O. 



5. The effective action of composite bosons plus 
qucisiparticles 

Finally I consider the case in which one needs to retain, in addition to 
composite bosons, also quasiparticles. To find the total effective action I 
first observe that the new coherent states can be rewritten 



la, b) = exp 



l-a{R*)iB^ R^a 



\a,b)) 



where I introduced the definition 

|a,6)) =exp(c'' R-^aj \b) . 



(40) 



(41) 



Notice that in the latter states there appear the original fermion operators, 
not the quasiparticle ones. 
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According to Eqs. (22,23) the action is given in terms of the ratios 

{aibi\a2b2)~^{aibi\exp{-TH)\a2b2) ■ (42) 
The matrix element of any operator O factorizes according to 

{ai,b^\0\a2,b2) = exp (-^alRl,B,{Rj,)i a*, - \a2{R*22)'-BlRl^a2 



x((ai,6i|0|a2,62)). 



(43) 



The first factor cancels out in the ratios (42) and the second factor is given 
by 

((ai,6i|exp(-T^F)|Q!2,62)) = exp(-T^F)|b2) exp(-/ - riJgp) . (44) 

The first matrix element yields the composite boson term Sc of the action 
reported in the previous Section, while the second one is the quasiparticle 
contribution. To shorten its expression I define the Grassmann variables 



w* = R 
w = R 

In terms of these variables 
1 



{R'^y~'a*-B'<R---a 
R-ia + B{R'^y^ a* 



(45) 



/ 



w'w — -w*Bw* — -wB^w 
2 2 



w* ih + Y,9KFi,R^FK \ w--Y,gKtr [(^B^ RF^ 

\ K J K 

X wFkw + {Fk RB) w*Fkw*] -J2gK ^w*f1.w* ^wFkw .(46) 

K 

Now I restrict myself to the case of a number of quasiparticles much smaller 
than 0. Then assuming Vt^K ~ ^t^x ~ 1 the noncanonical temporal terms 
in the action are of order Ct~^ 

i<(6t_i - Bt)w; + ^Wt-dBl - Bl,)w,_, = Oin--^) (47) 

because B = fl^^b ■ and b2 — bi = 0(1). To this approximation the 
temporal terms in the action become a| (ai — a2), which are canonical, and 
therefore the total Hamiltonian can be written in operator form 



H = Hb+H, 



IP 



(48) 
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where 



Hqp 



[. K K 

X wFkw + {FkRB) w^Fj^w^'^ -^9k ^w^Fr 




K 

In this expression the operators w'^ , w must be regarded as functions of the 
boson and quasiparticle operators ¥,b,a^,a and the colons mean normal 
order with respect to the latters. 

I do not have the space to discuss the properties of the Hamiltonian I 
have presented. I only make the following observations 

i) it respects term by term all fermion symmetries, in particular fermion 
number c;onservation, as it can be seen remembering that the bosonic op- 
erators w have fermion number 2 

ii) it contains quadratic quasiparticles terms of the type ce^ce^ and ad, 
similar to the "dangerous" terms of Bogoliubov. These terms however do 
not break fermion number conservation because they are accompanied by 
the operators b, respectively. Using the variational equations satisfied 
by the structure functions of the composites, Eq.(81) of Ref. 2, one can 
see that the coefficients of these terms vanish and the term involving the 
operator a^a takes the usual form. Therefore the difference with respect to 
the standard Bogoliubov transformation is that the quartic quasiparticle 
term involves interactions whith composite bosons. 
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